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We present the results of several parametrizations to two different ensemble of data on pp total 
cross sections (cf^ t ) at the highest center-of-mass energies (including cosmic-ray information). The 
results are statistically consistent with two distinct scenarios at high energies. From one ensemble 
the prediction for the LHC (y/s = 14 TeV) is a v t v ot = 113±5 mb and from the other, a v t v ot = 140±7 mb. 
From each parametrization, and making use of derivative analyticity relations (DAR), we determine 
■ p(s) (ratio between the forward real and imaginary parts of the elastic scattering amplitude). A 

' discussion on the optimization of the DAR in terms of a free parameter is also presented. In all 

cases good descriptions of the experimental data are obtained. 
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I. INTRODUCTION 



■ The total cross section, a to t , and the p parameter (ratio of the real to imaginary part of the forward elastic scattering 
^ | amplitude) are important quantities in the investigation of elastic hadron-hadron scattering at high energies jl]^] . They 

. are expressed in terms of the elastic scattering amplitude, F(s,t), by the formulas, 

vo 

O 1 <?tot{s) = 4arImF(a, t = 0), (1) 

in 

ReF(s,t = 0) 

P{S) = ImF(s,t = 0Y (2) 

where y/s is the center-of-mass energy and t the four-momentum transfer squared. 

For proton-proton (pp ) collisions, both quantities have been extracted from accelerator experiments and the results 
O || extend up to y/s ~ 63 GeV. On the other hand, of% t may also be inferred from cosmic-ray experiments at still higher 
energies, y/s ~ 10 TeV. However, these cosmic-ray results are model-dependent, since they are obtained from proton- 
air cross sections || and this has originated different results for of^ in the interval y/s = 5 — 25 TeV, which exhibit 
a reasonable degree of discrepancy. 

Several models present extrapolations at the cosmic-ray region and the observed discrepancies may be accommo- 
dated by different models {§-§]. As we shall recall, presently, it is difficult to decide which could be the "correct" 

■ result and it seems that the general trend is to expect the new data from the next accelerator experiments, the BNL 
Relativistic Heavy Ion Collider (RHIC), y/s~ ~ 500 GeV § and the CERN Large Hadron Collider (LHC), y/s ~ 14 
TeV ®. 

However, at this stage, we understand that a model independent analysis of the experimental information presently 
available, taking detailed account of the discrepancies and its consequences, may bring important insights on the 
subject. This is the central point we are interested in. 

In this communication, we first investigate two sets of distinct results for cr^ t , at cosmic-ray energies, in a model 
independent way that also takes into account of the experimental data at the accelerator region (y/s > 10 GeV). In each 
case we fit four different analytic functions to each ensemble of data, using the CERN-MINUIT routine |jll[ . Next we 
make use of the one-subtracted derivative analyticity relation (DAR) in order to obtain analytic parametrizations for 
the p parameter as function of the energy, from both ensembles and for all the parametrizations. In this investigation 
we stress that the general expression of the DAR has a free parameter and we present a study on the practical 
applicability of this parameter in the fits to the p data. In all cases good descriptions of the experimental data are 
obtained. 

The paper is organized as follows. In Sec. II we review the experimental information on a^ t from accelerator and 
cosmic-ray experiments and recall the main steps connecting integral and derivative analyticity relations. The fits 
concerning af^ t and the results for p(s) are presented in Sec. III. In Sec. IV we discuss in some detail all the results 
obtained and present comparisons with some model predictions. The conclusions and some final remarks are the 
content of Sec. V. 
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II. EXPERIMENTAL INFORMATION AND ANALYTICAL APPROACH 



Here we first review the experimental information presently available on af^ t at the highest energies, stressing the 
puzzles involved at the cosmic-ray region. We also recall some essential formulas connecting integral and derivative 
analyticity relations and the corresponding derivative relation between u\^ t and p (the analytical approach) . 



A. Experimental information on pp total cross sections 

As mentioned before, in order to analyze the experimental information presently available on <Jtot{s) for pp inter- 
action at energies beyond ~ 10 GeV (high energy region), we must distinguish between accelerator and cosmic-ray 
information. In the later case we follow a discussion first presented in Ref. [Q. 

Concerning accelerators, data on of^ from three experiments at the CERN Intersecting Storage Ring (ISR), between 
23.5 GeV and 62.3 GeV, were critically analyzed by Amaldi and Schubert and we shall consider here the final mean 
values from their analysis fT^| . We also include the results at 13.8 GeV and 19.4 GeV, obtained in Fermilab fl2|| . 

Although other accelerator data exist in this region, this set is suitable for the analyzes we are interested in. From 
one hand, these data represent the correct trend of pp collisions in the region 10 — 60 GeV and allow a statistical 
approach that includes different cosmic-ray information (see what follows). On the other hand, as will be discussed in 
Sec. IV, they are adequate for comparison with a model that predicts a faster rising of af% t than generally expected. 

Differently from accelerator data, cosmic-ray experimental information on af^ t comes from proton-air cross sections. 
We recall that antiprotons are not expected to have any significant role in these interactions and therefore the bulk 
of cosmic ray information on hadron-nucleus cross sections does not concern antiprotons. 

Now, since what is extracted in these experiments is the proton-air cross section, the determination of the pp cross 
section depends on nuclear model assumptions and this has originated some puzzles, as reviewed in what follows. 

The first information on a^ t from cosmic ray experiments, at y/s = 30 TeV, was reported by Baltrusaitis et al. in 
1984 fll3|] . Extensive air showers recorded by Fly's Eye detector allowed the determination of the proton-air inelastic 
cross section. Based on Glauber theory, assuming a Gaussian profile for the nucleus and the proportionality between 
total cross section and the slope parameter (geometrical scaling), the authors inferred Jl3|] 

a tot = 120 ± 15 mb at y/s = 30 TeV. (3) 

In 1987, based also on data from Fly's Eye experiment, Gaisser, Sukhatme and Yodh (GSY) introduced the limit 
&tot ^ 130 mb at y^s ~ 30 TeV and, taking account of various process in the Glauber theory and additional 
assumptions, estimated Q 



a 



fP t = 175+^0 mb at ^ = 40 TeV. (4) 

Afterwards, extensive air shower data collected in the Akeno observatory allowed new estimates in the range 5 — 25 
TeV, reported by Honda et al. |u| , which is in agreement with the result reported by Baltrusaitis et al. and in 
disagreement with the GSY result. In particular, extrapolations by Honda et al. indicated [ fi"5| 

aZ = 133 ± 10 mb at Vs = 40 TeV, (5) 

which is in agreement with the result reported by Baltrusaitis et al. and in disagreement with the GSY result. In 
the same year, Nikolaev claimed that the Akeno results have been underestimated by about 30 mb |ll| and, if we 
take this correction into account, the data in the interval 5 — 25 TeV show agreement with the early GSY result. An 
important point is the fact that the analysis by Nikolaev seems correct, has never been criticized and the same is true 
for the GSY limit and result. 

All these cosmic-ray informations are shown in Figure 1, together with the accelerator data at lower energies. We 
stress that these experimental information concern only pp and not pp interactions. From this Figure we clearly see 
the discrepancies between Honda/Baltrusaitis and GSY/Nikolaev. 

For these reasons, in this paper we investigate the behavior of the total cross section taking account of the dis- 
crepancies that characterize the cosmic ray information. In order to do that, we consider two ensembles of data and 
experimental information, with the following notation: 

Ensemble A: accelerator data + Honda et al. + Baltrusaitis et al. 

Ensemble B: accelerator data + Nikolaev + GSY 

For ensemble A the references are @|]|[n§ and for ensemble B, In Sec. Ill we make use of different 

parametrizations to fit the data from each ensemble. 
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FIG. 1. Experimental information on pp total cross sections: accelerator data in the interval 13.8 < y^s < 62.5 GeV and 
cosmic-ray results in the interval 6.3 TeV < ^/s < 40 TeV. 



B. Derivative analyticity relations and the p parameter 



Both cr tot (s) and p(s) play a central role in the investigation of high-energy hadron scattering. Due to the connections 
between forward real and imaginary parts of the scattering amplitude, Eqs. (1) and (2), the analyticity (dispersion) 
relations constitute a suitable, model independent, approach for a simultaneous study of these quantities. 

Although integral dispersion relations have been widely used in the study of hadronic scattering, in general, the 
analytical and/or numerical integrations are not an easy task. However, at sufficiently high energies, the smooth 
increase of of^ (Fig. 1) allows to connect the integral form with a derivative one, which is easier to handle. The so 
called derivative analyticity relations (DAR) were introduced in the seventies jl7| and since then have been critically 
investigated jl8|,[l9| . Recently, a recursive approach was developed, as well as generalizations to an arbitrary number 
of subtractions, for both cross even and odd amplitudes, near the forward direction |2(}]. As in the integral case, the 
convergences may be controlled by subtractions and specific formulas are associated with cross even and odd functions 
(scattering amplitudes in the case of particle-particle and anti particle-particle interactions) ^ . 

At this point we stress once more that we are only interested in pp scattering (where discrepancies happen in atot(s)) 
and not pp. Besides, our investigation concerns the highest energies, characterized by the smooth increase of trf^ , 
nearly as a power on Ins. For these reasons in our analyzes we will consider only an even amplitude, as the leading 
contribution, and only one subtraction. 

A detailed deduction on how to obtain DAR from integral relations may be found in references jlT] . Here we 
only review the main steps concerning our case of interest, namely, one subtraction and an even amplitude. We begin 
with the well known once-subtracted integral dispersion relation (even amplitude) in the forward direction (t = 0) 

§13 

Ref+(s)=K+—I, (6) 
where K = Ref+(0) is the subtraction constant and 
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f +00 -1 
I = P I da? 



s'(s' 2 - s 2 ) 



-ImU(s'). 



(J) 



Following Bronzan, Kane and Sukhatme |17fl , and also pOj| . wc consider a real parameter a so that after multiply 
and divide by s a and integrating by parts we obtain 



I = — In 

2ss' 



— I dsV a-1 ]n 
2s 



(S) 



Taking account of the high-energy region (« > so ~ m2 ~ 1 GeV 2 ) and performing a change of variable s' = , s = e*, 
the last equation may be put in the form 



I= hi°° d ^' s ' a ' 1 lncoth \^ " ei { a ~ 1 + ^) Im M s 'V s " 



(9) 



Expanding im/ + (s')/ s' a in powers of £' — £, after some manipulation and taking account of the high energy limit 
(sq — > 0, that is, In s ~~ * —00) we obtain 



Ref + (s) =K + S «jr (Jm/ + M/««) § 



n=0 



where /„ represents the integral in the variable 

J„ = i | + J ^'e^- 1 )^'-^ lncoth ||r - £| (a - 1 + (£' - 0' 



(10) 



(11) 



We have assumed that the series may be integrated term by term. Denoting £' — £ = y and integrating by parts, 
this equation may be put in the form 



/„ = - lncoth e^Vll* + - / dy e ——y n . 
ir 2 00 7r„/-oo smhy 

A central point here is that the first term on the right hand side converges to zero only for 

< a < 2. 



(12) 



(13) 



In this case we only have the second term, which can be expressed as a recursive relation in terms of the parameter 
a |0|: 



da( n ) 



and integration in the complex plane gives 



I a = tan Q(a - 1)) . 



(14) 



(15) 



With this, we obtain an expression connecting the real part of an even amplitude with the derivatives of the 
imaginary part at the same energy, namely the DAR: 



Ref + (s) K 
= h tan 

s a s a 



-\a-l 



d In s 



The leading term in the tangent series reads 

/TT , 

tan — [a — 1 
\2 y 



Imf + ( 3 ) n 2f n 
s a 2 



Imf+{s) 



d ( Imf+(s) 



d\n s 



(16) 



(17) 
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and making use of the normalization (fc 2 ~ s), 

/+(«) _ 
s 

we obtain the general result for a forward amplitude: 

K 



F(s,t = 0), (18) 



ReF(s, 0) = h tan (- (a - 1)) ImF(s, 0) + 

„ , 7r o/7r, s \ d ( ImF(s,0)\ 



d In s \ s Q 1 

At last, from equation (1), we obtain the general relation connecting p and o~tot'- 

, . 4irK /7T. a 7r s Q_1 o/7r, d /c tot (s)\ ,„„. 

P S = rr + tan - a - 1 + --sec 2 - a - 1 ) — . 20 

scrtot(s) V2 / 2 <T tot {s) \2 ) dins \ s / 

The standard form referred to in the literature and applied to hadronic scattering corresponds to the particular 
choice a — 1 (but usually without the subtraction constant) flp| p2||2^| , which we shall refer as the conventional form 
of the DAR: 

. , AirK 1 7T d , , .. , „ 

scrtot(s) cr tot (s) 2 dins 

In the next section we use both the conventional and the general relations, in order to determine p(s) from different 
parametrizations for a^ t (s) and ensembles of data. 

III. FITS AND RESULTS 

In this section we first present the fits to the total cross sections for both ensembles defined in Sec. II and then the 
predictions for p(s) obtained by means of DAR. The discussion on all the obtained results is the content of Sec. IV. 



A. Fits to total cross section 

For each ensemble defined in Sec. II, we fit the data with some standard and suitable parametrizations for the total 
cross sections: 

fit I: a tot = A + Blns + C{liis) 2 (22) 
fit 2: atot = A + B\ns + C(\ns) D (23) 
fit 3 : a tot = A + B In s + C(ln s) 2 + i?s~ 1/2 (24) 

fit 4: cr fot = yl + Slns + C*(lns) D + i?s- 1/2 (25) 

where A, B, C, D and R are free parameters. 

The choice for these parametrizations was based on the following considerations. Firstly, since Figure 1 is a linear-log 
plot, we see that at the highest energies the data suggest an increase of af^ t as a polynomial on Ins. The Froissart- 
Martin bound states that the fastest rate permissible asymptotically for the rising of af^ t is ln 2 s |2^]. However, it 
has been shown by the UA4/2 Collaboration that fits to pp and pp accelerator data indicate the power 2. 251° ^, 
a result referred to as a "qualitative saturation of the Froissart-Martin bound" |2j|. For these reasons, we consider 
polynomial functions on Ins with two possibilities for the power factor: the value 2, according to the asymptotical 
bound and as a free parameter to be determined by the fits. Finally, the power function on s, Eqs. (24) and (25), 
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represents the usual way to take account of data at lower values of the energy (5 ~ 20 GeV) and can be associated 
with Regge phenomenology. 

The fits were performed by using the CERN-MINUIT routine and the results for ensembles A and B are shown 
in Figs. 2 and 3, respectively. The corresponding values of the free parameters and the x 2 per degree of freedom are 
displayed in Table 1 and 2 for ensembles A and B, respectively. We discuss all these results in Sec. IV. 

Table 1. Values of the parameters in Equations (22) — (25) and the x 2 per degree of freedom in each fit to 

Ensemble A. 



fit: 


1 


2 


3 


4 


A 


45.78 


42.64 


38.74 


33.67 


B 


-3.315 


-1.875 


-1.868 


-0.3201 


C 


0.3654 


0.1220 


0.2894 


0.09465 


D 


2 (fixed) 


2.312 


2 (fixed) 


2.288 


R 


- 


- 


21.52 


31.40 


x 2 


11.7 


15.4 


9.8 


9.9 


d.o.f. 


11 


10 


10 


9 


X 2 /d.o.f. 


1.06 


1.54 


0.98 


1.10 


Table 2. 


Values of the parameters 


in Equations (22) — (25) and the x 2 


per degree of freedom 


in each fit to 






Ensemble B. 






fit: 


1 


2 


3 


4 


A 


49.27 


43.54 


62.48 


38.34 


B 


-4.435 


-1.846 


-7.221 


-0.8162 


C 


0.4534 


0.05958 


0.6050 


0.03461 


D 


2 (fixed) 


2.624 


2 (fixed) 


2.755 


R 






-38.40 


14.80 


x 2 


15.0 


10.7 


10.4 


11.4 


d.o.f. 


11 


10 


10 


9 


X 2 /d.o.f. 


1.36 


1.07 


1.04 


1.27 




FIG. 2. Fits to pp total cross section data from ensemble A through Eqs. (22) - (25): fit I (solid), fit 2 (dotted), fit 3 
(dot-dashed ) and fit 4 (dashed). 




FIG. 3. Fits to pp cross section from Ensemble B. Same legend as Fig. 2. 
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B. Predictions for the p parameter 



In this section we determine the p(s) behavior, making use of the DAR and all the parametrizations for <rf^ t obtained 
in the last section. As commented before we are interested in the practical applicability of both forms of the DAR, 
the general and the conventional ones, Eq. (20) and Eq. (21), respectively. 

In principle, the general expression (20) has two "free parameters" , namely the subtraction constant K and the 
parameter a. Firstly, as a simple exercise to see the effect of the DAR, let us consider one of the fits to of^t > f° r 
example fit 4, for both ensembles A and B, and calculate p{s) through the conventional expression and without 
subtraction constant, which means taking a = 1 and K = in Eq. (20), or K — in Eq. (21) (this corresponds to the 
formula usually referred to in the literature Jl|,^| ) . The results are displayed in Fig. 4 together with the experimental 
data [^6|. We see that although the predictions from both ensembles are similar at the ISR energy region, both 
disagree with the data. 

In what follows we investigate the influence of the above two free parameters in the description of these data. In 
order to do that, we first treat the conventional formula with the subtraction constant and then the general formula, 
without subtraction constant and the factor a as free parameter. 

• Conventional derivative dispersion relation 

With the parametrizations (22-25) for o~\^ t we obtain analytical expressions for p(s) by using the conventional form 
of the DAR, Eq. (21). For each input parametrization the subtraction constant if is a free parameter, determined 
by fit, through the CERN-MINUIT routine, to the experimental p data 

The results from each parametrization and for both ensembles are shown in Figures 5 and 6 together with the 
experimental data. The values of the subtraction constant and statistical information about the fits are displayed in 
Table 3. 



0.20 - 




FIG. 4. Results for p(s), using fit 4 from both ensembles A (dashed) and B (solid) and DAR with a — 1 and K — in Eq. 
(20) and experimental data pif . 
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Table 3. Values of the subtraction constant from fits to p(s) data and the % 2 for 6 degree of freedom. 
Calculation performed thorough Eq. (21) using the four fits to of^ t from ensembles A and B. 



fit-of„ t - Ensemble A: 


1 


2 


3 


4 


K 


-134.1 


-135.6 


-130.3 


-129.6 


X 2 /d-o.f. 


0.82 


0.63 


0.97 


0.63 


fit-<r££ - Ensemble B: 


1 


2 


3 


4 


K 


-130.0 


-134.4 


-137.1 


-132.0 


X 2 /d.o.f. 


2.0 


1.48 


1.85 


1.44 





0.20 - 




FIG. 5. Results for p(s), through the conventional DAR, Eq. (21), from fits to af£ t from ensemble A together with the 
experimental data. Same legend as Fig. 2. 
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0.20 - 




FIG. 6. Results for p(s), through the conventional DAR, Eq. (16), fits to af£ t from ensemble B, together with the experimental 
data. Same legend as Fig. 2. 



General derivative analyticity relation 



As shown in some detail in Sec. II. B, the general result for the DAR depends on the free parameter a. It comes 
from the integration by parts of Eq. (7) and it is necessary in order to allow a finite derivative form associated with 
the integral form. In this sense it seems to play the role of a regularization factor. Although in nearly all practical 
uses of the DAR the value a = 1 is assumed, we have shown that its value is constrained to the interval < a < 2 
(see also references |l^,|l9|). In particular, in the context of a multiple diffraction model, it was recently shown that 
the description of the experimental data on pp elastic scattering may be improved by taking a as a free parameter 
f27f . This early result inspired us to make use of the general relation (20) and to investigate the possible effect of a as 
a free parameter (some preliminary results were already presented in Ref. p8[|). To this end we will not take account 
of the subtraction constant in the general formula for the DAR, namely K = in Eq. (20), so that we explicitly have: 



/ o(s)=tan(|(a-l)) -+ 
(=(-!) 



— sec 
2 



1 da t Qt(s) 
ctot(s) dins 



1 - a 



(26) 



With a fixed parametrization for a\^ t we can fit the experimental data on p by letting a to be a free parameter in 
the above equation, once more by using the CERN-MINUIT. The results from all the parametrizations for af^ t with 
both ensembles A and B are shown in Figures 7 an 8, respectively. The corresponding values for a and statistical 
information about the fits are displayed in Table 4. 



Table 4. Values of the parameter a from fits to p(s) data and the \ 2 f° r 6 degree of freedom. Calculation 
performed thorough Eq. (20) with K = and using the four fits to of^ from ensembles A and B. 





fit-of p t - Ensemble A: 


1 


2 


3 


4 


Q 

X 2 /d-o.f. 


1.229 
2.61 


1.223 
3.17 


1.231 
2.17 


1.230 
2.20 


frt-a™ t - Ensemble B: 


1 


2 


3 


4 


a 

X 2 /d-o.f. 


1.244 
1.48 


1.240 
1.94 


1.244 
2.00 


1.239 
1.80 
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0.20 : 
0.15 - 




FIG. 7. Results for p(s), through the general expression of the DAR, Eq. (20), with a as free parameter, K = and fits to 
afL from ensemble A. 



0.20 - 




FIG. 8. Results for p{s), through the general expression of the DAR, Eq. (20), with a as free parameter, K = and fits to 
a^ t from ensemble B. 
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IV. DISCUSSION 



The parametrizations for of^ t are the usual ones, but the different ensembles, suggested by cosmic-ray results, 
introduce novel behaviors in the asymptotic region, as clearly shown in Figs. 2 and 3. For example, by calculating 
the average value and the standard deviation from the four parametrizations we can estimate cf^ (14 TeV) = 113 ±5 
mb for ensemble A and a^ t (14 TeV) = 140 ± 7 mb for ensemble B. At lower values of the energy, namely ^fs 
in the region 10 — 100 GeV, there is no significant distinction between the four parametrizations. In general, 
model predictions are in agreement with the above result from ensemble A J^-Q], including the fit by the UA4/2 
Collaboration j25|. To our knowledge, the only exception (model) that presents agreement with the result from 
ensemble B is that of Ref. || . This multiple diffraction model is based on analyzes of pp elastic scattering in the 
interval 13.8 GeV < y/s < 62.5 GeV, the same set we used here at the accelerator region. Extrapolation to higher 
energies predicted <7t t(16 TeV) = 147 mb, without estimated errors. Recently Perez-Peraza et al. improved this 
model predictions determining confident error bands through a forecasting regression analysis |p9f . Reading from Fig. 
2 of this reference we can infer cr tot (16 TeV) ~ 147 ± 37 mb. Despite of the large error band ( ~ 25%), even in this 
case it is clear from the quoted figure that the results favors ensemble B. From a "statistical point of view" we may 
say that the published results from models and fits show agreement with ensemble A. However, it should be stressed 
that the cosmic-ray estimations in the ensemble A were obtained under the hypothesis of the geometrical scaling 
|3j, which is violated even at \fs ~ 500 GeV. Moreover, we should remember results from cosmic-ray experiments 
which indicate the possibility of new phenomena in pp collisions at center-of-mass energies beyond 500 GeV ||cj . As a 
direct consequence we should expect new open channels and therefore a faster rising of the pp total cross section than 
expected in the extrapolations from accelerator data. This seems to be well accommodated by the parametrizations 
with ensemble B. 

Concerning the determination of p(s) through DAR, we first showed that the conventional expression without 
subtraction constant (a — 1 and K = 0) does not reproduce the experimental p data (Fig. 4). We recall that this 
expression has been referred to in the literature 
Taking account of the subtraction constant K, t 



l]||] and also used in the context of phenomcnological models p2| . p3[ 
re data are well described, as shown in Figs. 5 and 6. From Table 3, 



\K\ ~ 130 - 137, (27) 

depending on the fit and ensemble for af^ t . As expected different ensembles correspond to distinct behaviors at the 
highest and asymptotic energies. As in the case of a\^ t we can estimate an average value and standard deviation from 
the four results in Figs. 5 and 6: p(14 TeV) = 0.142±0.010 with of G p t from ensemble A and p(U TeV) = 0.173±0.013 
from ensemble B. In this case we can say that the results with ensemble A present the best agreement with the 
experimental data (Table 3 and Figs. 5 and 6). We also tested the general expression for the DAR by taking K = 
and letting a as a free fit parameter. The results presented in Figs. 7 and 8 show that the data are also satisfactorily 
described, specially in the case of af% t from ensemble B (compare with Fig. 4, the case of fixed a = 1). From Table 4, 

a - 1.22 - 1.24, (28) 

depending also on the fit and ensemble for af^ t . In all these cases the condition (13) is verified. 



V. CONCLUSIONS AND FINAL REMARKS 



In this communication we investigated two ensembles of experimental information on pp total cross sections and 
used four different parametrizations to fit the data, as function of the energy. In each case we obtained predictions 
for p(s) making use of both the conventional and the general expressions for the derivative analyticity relations. 

Our first main conclusion is that experimental information presently available on pp total cross sections indicates two 
possible different scenarios for the hadronic interactions at the highest energies. The fast rise of a\^ t from the analysis 
of ensemble B is corroborated by the multiple diffraction model of Ref. g and also indication of new phenomena from 
emulsion chamber experiments p0| ]. Although in this work we only point out for this possibility, new information 
coming from RHIC and LHC shall certaintly clarify this subject. In this sense, our results may be viewed as a kind 
of warning against some possible precipitated assumptions, namely, that extrapolations from accelerator data which 
show agreement with the cosmic-ray estimations in ensemble A could be the final answer to the question. 

A second novel result from this model independent analysis was to show the practical applicability of both the 
conventional and the general expressions for the derivative analyticity relations at sufficiently high energies. In the 
conventional case, a = 1, letting the constant K as free parameter, the description of the experimental p data is quite 
good. On the other hand, taking K = and letting a as free fit parameter, the results may be considered satisfactory. 
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In both cases the predictions at the highest energies are practically the same (compare Figs. 5 and 7 and also Figs. 
6 and and 8) and the differences at this energy region come obviously from the different ensembles for <j\^ t . Our 
approach was to consider the two possibilities separately (A' ^ or a ^ 1), so that we could infer the intervals of 
possible variations, Eqs. (27) and (28). Simultaneous analysis with both possibilities shall improve the description of 
p(s) and this is our second main conclusion. 

In this communication we treated only pp interactions, since the cosmic-ray informations concern only this case. 
We are presently investigating the inclusion of antiproton-proton data in the analysis through adequate considerations 
on crossing symmetry. 

We are grateful to Fapesp for financial support and J. Montanha for reading the manuscript. M. J.M. is also thankful 
to N.N. Nikolaev for discussions, V. Ezhela for correspondence and CNPq for financial support. 
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